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Abstract 

Many aspects of high-temperature gauge theories, such as the electroweak 
baryon number violation rate, color conductivity, and the hard gluon damp- 
ing rate, have previously been understood only at leading logarithmic order 
(that is, neglecting effects suppressed only by an inverse logarithm of the 
gauge coupling). We discuss how to systematically go beyond leading loga- 
rithmic order in the analysis of physical quantities. Specifically, we extend to 
next-to-leading-log order (NLLO) the simple leading-log effective theory due 
to Bodeker that describes non-perturbative color physics in hot non-Abelian 
plasmas. A suitable scaling analysis is used to show that no new operators 
enter the effective theory at next-to-leading-log order. However, a NLLO 
calculation of the color conductivity is required, and we report the resulting 
value. Our NLLO result for the color conductivity can be trivially combined 
with previous numerical work by G. Moore to yield a NLLO result for the hot 
electroweak baryon number violation rate. 



I. INTRODUCTION 



The near equilibrium dynamics of hot, weakly-coupled non-Abelian plasmas, such as high 
temperature QCD or electroweak theory, involves a surprisingly rich hierarchy of spatial and 
temporal scales. Of particular interest in recent years is the analysis of non-perturbative 
processes occurring in (or near) equilibrium in such plasmas. A primary application is to the 
rate of baryon number violation in hot electroweak theory, which quite possibly is responsible 
for the preponderance of matter over anti-matter in our universe [1]. Baryon number non- 
conservation is possible due to the electroweak anomaly in the baryon (and lepton) number 
current, and proceeds through non-perturbatively large thermal fluctuations in electroweak 
gauge fields. A key goal is to understand how to calculate the dynamics of such fluctuations. 

As will be reviewed momentarily, non-perturbative gauge dynamics in hot non-Abelian 
plasmas may be described (at leading logarithmic order) by a remarkably simple effective 



theory discovered by Bodeker [2], whose only input parameter is the value of the color 1 
analog of electrical conductivity. The color conductivity in turn depends on what is known 
as the hard gauge boson damping rate, which is a measure of the rate at which collisions 
randomize a charge carrier's color charge. Up to now, however, these quantities have only 
been known at leading order in the logarithm of the gauge coupling; that is, neglecting 
effects suppressed only by 1/ ln(l/g), where g is the gauge coupling. Consequently, numerical 
simulations of the high temperature topological transition (or baryon number violation) rate 
based on Bodeker's effective theory have only been valid to leading order in logarithms [3]. 
For practical purposes such leading log results are, by themselves, not very useful — there is 
a huge difference between, say, \n(lQir/g) and ln(l/10g) for any realistic value of the gauge 
coupling. But this difference is of sub-leading order in ln(l/g). In this paper, we report a 
calculation of the color conductivity to next-to-leading-log order (NLLO), and in the process 
explain what this means in gauge-invariant language. Our approach involves systematically 
constructing a sequence of effective theories which reach from short to long distance scales, 
and performing perturbative calculations which enable one to match the physics of each 
effective theory with its predecessor to the order desired. When combined with Moore's 
numerical result [3] for the constant of proportionality between the topological transition 
rate and the color conductivity, our improved value for the conductivity immediately yields 
a NLLO result for the topological transition rate. 

In this paper, "hot" plasma means hot enough (1) to be ultra-relativistic, (2) to ignore 
chemical potentials, (3) for non-Abelian gauge couplings to be small, and (4) to be in the 
high temperature symmetric phase. 2 For convenient reference, the critical distance and time 
scales associated with color dynamics in such plasmas (to be reviewed below) are collected 
in table I, and the notation used in this paper is summarized in table II. 



0(T) typical momentum of excitations in the plasma. 

m = O(gT) inverse Debye screening length. 

7 = 0[g 2 T \n(g~ 1 )] hard gauge boson damping rate. 

0(g 2 T) inverse spatial scale of non-perturbative gauge fluctuations. 

0[g 4 T In (g^ 1 )] inverse time scale of non-perturbative gauge fluctuations. 



TABLE I. Important (inverse) distance and time scales. 



1 We use "color" as a generic term for some non-Abelian gauge charge, not as something specific to QCD. 

2 We use the term "symmetric phase" loosely since, depending on the details of the Higgs sector, there 
need not be any sharp transition between the symmetric and "symmetry-broken" phases of the theory [4] . 
A sharp transition is required for viable electroweak baryogenesis. Our analysis is applicable whenever the 
infrared dynamics of the Higgs field is irrelevant at lengths of 0(1/ g 2 T), which is the case either (a) far 
above the electroweak phase transition or "crossover," or (b) in the symmetric phase at the transition in 
cases where there is a sufficiently strong first-order transition. 
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t> M = (l,v); v a spatial unit vector. 

A = A(x, t), the spatial non-Abelian gauge field. 

D = V + gA a T a , the gauge covariant derivative. 

W = W(xl, v, t), the adjoint color distribution of hard excitations. 

C = C( x ; t) and £ = £(x, v, t) are Gaussian white noise. 

((■■■)) denotes averaging over noise. 

(•••} = (•• -) v denotes averaging over the direction v. 

is a 5-function on the two-sphere normalized so that (<5 52 (v— v')) v / = 1. 
5CW = (6C(v, v') W(v')) v /, the linearized collision operator applied to W. 
Ca is the adjoint Casimir of the gauge group [N for SU(iV)]. 
d = 3 — e with e — > 0, the number of spatial dimensions. 

TABLE II. Summary of notation. 
II. REVIEW 

It is well known that non-perturbative gauge field fluctuations in hot non-Abelian plasmas 
are essentially magnetic and characterized by the distance scale 3 R ~ (g 2 T)~ 1 . The corre- 
sponding time scale has only been understood more recently, and can be arrived at by the fol- 
lowing simple physical argument. Start with (non-Abelian) Ampere's Law, D x B = D t E+j. 
Now recall that plasmas are conductors, and so write j = <rE, where a is the (color) conduc- 
tivity. If we assume the time dependence will be slow and treat the covariant time derivative 
D t E as negligible compared to erE (an approximation that can be verified a posteriori) , then 
the result is 

DxB = (jE. (1) 
Extracting the time scale from this equation is clearest if one specializes to A = gauge: 

D x D x A = a— A . (2) 

dt 

In terms of characteristic scales, this is just 

^A-a^A, (3) 

and so the characteristic time t ~ aR 2 ~ a/g A T 2 . The color conductivity o was originally an- 
alyzed by Selikhov and Gyulassy [7] and is order a ~ T j ln(l/g), yielding t ~ [(^Tln^/g)] -1 . 
The rate per unit volume T for topological transitions in a hot, non-Abelian plasma must 
therefore have the characteristic size T ~ l/(R 3 t) ~ a 5 T 4 ln(l/a), where as usual a = g 2 /47r. 
The factor of a 5 in this rate was originally noted by ourselves and Son [8], and the presence 
of the logarithm by Bodeker [2]. The explanation in terms of color conductivity is given in 
detail in Ref. [5]. 



3 See the introduction of Ref. [5] for a simple heuristic argument, or Ref. [6] for more formal ones. 
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The above order of magnitude arguments can be made quantitative. One can nu- 
merically integrate, with Boltzmann-weighted initial conditions, the classical equation (2) 
that we derived from Ampere's Law, and measure the frequency of topological transitions 
[3]. A classical, as opposed to quantum, treatment is adequate because the low momen- 
tum modes (k <C T) of the gauge field in a hot plasma have large occupation numbers 
[nk ~ (e fc / T — ~ T/k ^> 1] and so, by the correspondence principle, behave classically 
[9]. Corrections to the classical approximation are suppressed by powers of k/T. In our 
case, we are interested in modes with k ~ i? -1 ~ g 2 T, which is indeed small compared to T 
for weak coupling. 

However, Eq. (2) first requires some refinement. It is a purely dissipative equation. An 
effective theory for the equilibrium dynamics of the system must therefore incorporate a 
source of thermal noise to maintain thermal equilibrium. Ultimately, this noise is (like the 
dissipation) due to interactions of the modes of interest with other degrees of freedom of the 
system, specifically typical charge carriers, which have momentum of order T. As discussed 
in Ref. [5], one may argue from general principles that, at the distance and time scales of 
interest, this noise can be taken to be simple Gaussian white noise, appropriately normalized 
to produce the desired temperature T. The resulting effective theory is 4 

DxB + ( = (tE, (4a) 

where £ is Gaussian noise with correlation 

((0(*,x)0(t',x')» = 2aTM(t-t'M (3) (x-x / ). (4b) 

The effective theory (4) is also known as stochastic 3-dimensional gauge theory, and has a 
long history of theoretical study motivated by formal questions divorced from any connection 
with high temperature dynamics [14]. A key property of this theory is the fact that it is UV 
finite [5]. 

The final ingredient needed to make the effective theory quantitative, is the precise value 
of the color conductivity a. This has previously been known only at leading log order [7,2,5]: 

2 

m 

^37' (5) 

where ~ denotes equality up to relative corrections suppressed by [ln(l/(?)] _1 . Here, m is 
the leading-order Debye screening mass, which is well known, 5 is O(gT), and completely 



4 This equation certainly makes sense in Aq = gauge, and is also correct and unambiguous in general 
flow gauges of the form A = R[A], where R[A] depends on A(x, t) only instantaneously and so does not 
involve time derivatives of A. See Ref. [10] for a discussion of flow gauges, and Ref. [11] for a proof that 
the equation (4) may be applied in any gauge of this class. There are subtleties in directly interpreting the 
Langevin equation (4) in other gauge choices, such as Landau gauge. For various historical reasons, our 
approach for doing actual calculations is to use a path integral representation of these Langevin equations, 
for which it is straightforward to fix any gauge desired using the usual Faddeev-Popov procedure. (See Ref. 
[12] for details.) We find it convenient, in particular, to use Coulomb gauge. But Coulomb gauge might also 
be implemented as the A — > oo limit of the flow gauge Aq = AV • A, for which there should be no problem 
working directly with the Langevin formulation (4) instead of the associated path integral representation. 

5 For SU(iV) gauge theory with N s scalars and JVf Dirac fermions in the fundamental representation, the 
Debye mass m 2 = ^(2N+N S + N{) g 2 T 2 at leading order in g. For hot electroweak theory with three fermion 
families and a single Higgs doublet, m 2 = ^ g 2 T 2 . 
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encapsulates the dependence on the matter field content of the underlying theory. 7 is the 
hard gauge boson damping rate [15] and, to leading-log order, is given by 




(6) 

The leading-log effective theory represented by (4) and (5) was first derived by Bodeker 
[2]. (For subsequent derivations, see Refs. [5, 3, 16-18]. 6 ) Using this theory, Moore [3] has 
numerically simulated the topological transition rate for electroweak theory, obtaining 

r « (10.8 ± 0.7) (fj)(«T) 5 . (7) 



III. THE SEQUENCE OF EFFECTIVE THEORIES 

The series of effective theories we will need to describe color dynamics at low frequency 
uo and momentum k [ultimately uo ~ t^ 1 ~ [g 4 T 'ln(l / g)]^ 1 and k ~ R^ 1 ~ (g 2 T)^ 1 } are as 
follows. 

Theory 1: uo,k<^T 

(D t + v • D) W - v • E = , (8a) 
D v F' a/ =f = m 2 (v^W) . (8b) 



Theory 2: 00 <C k <C m 

v-DFU-vE = -5CW + Z, (9a) 

(W) = 0, (9b) 

D x B = j = m 2 (vW) , (9c) 
2T - 

««» = — SC. (9d) 



Theory 3: to k ^ 7 

(jE = DxB + C, (10a) 
((CO) =2aT. (10b) 



6 The derivation in Ref. [16] treats the color charge of the charge carriers as classical, using Wong equations. 
This approach has been used before in analyzing hot gauge theories, and it is worth pointing out one way 
to understand, after the fact, why it works. In the case at hand, note that the result (5) for a does not 
depend on the color representations of the charge carriers in any way except implicitly in the value of m 2 . 
At leading order, the contribution to m of any charge carrier species is proportional to the adjoint Casimir 
Cr for the color representation of that species, and there is no other representation dependence. If one 
extracts the proportionality constant for representations with arbitrarily large values of Cr, one will then 
also automatically have the right proportionality constant for small representations. But the color of large 
representations can be treated classically. So treating the color of the charge carriers classically and then 
setting the square of the classical charge to the quantum value Q 2 = Cr at the end of the day works because 
of the simplicity of the color structure for the quantities of interest. 
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All theories have been written in forms that are local in space and time. All fields should 
be interpreted as classical. 

Theory 1. The first effective theory is a now-standard formulation of "hard-thermal- 
loop" dynamics [19] which is valid for iv,k <C T, up to corrections suppressed by powers 
of g. These equations amount to linearized, collisionless, non-Abelian, Boltzmann-Vlasov 
kinetic theory. One has conceptually split the degrees of freedom in the underlying quantum 
field theory into those associated with quanta that have large momenta p ~ T and those with 
low momenta k < T. The high momentum (or "hard") quanta are treated collectively by a 
linearized Boltzmann equation; the (bosonic) low momenta modes, referred to as "soft" (and 
sometimes also semi-hard) modes are, because of the high occupation numbers mentioned 
earlier, treated as comprising a classical field. The first equation (8a) is a collisionless 
linearized Boltzmann equation for the propagation of the hard quanta in the presence of 
the soft gauge fields. W(~x.,v,t) represents the color distribution of hard particles in space, 
velocity, and time, where the velocity v is a unit vector because all hard particles are 
ultra-relativistic, and where (• • •) in Eqs. (8) and (9) denotes averaging over the direction 
of v. Although the hard particles have individual momenta p ~ T associated with very 
short wavelengths, the collective distribution of hard particles, described by W, can have 
slow spatial variation (in other words, the density of hard particles may vary slowly over a 
large region of space), which is why W appears in the effective theory describing u,k < T 
dynamics. 

The distribution W appearing in the above effective theories only encodes information 
about the color structure of collective fluctuations, and W lives in the adjoint representation 
of the gauge group. 7 (See Refs. [5,2,20] for details.) A full description of all physics for 
iv, k <C T would also contain additional distribution functions describing fluctuations that 
do not contribute to the color current. These include color-neutral fluctuations which are 
responsible for hydrodynamic phenomena such as viscosity and sound, but which do not 
couple to the long-distance color dynamics except at higher order in the coupling g. We 
will therefore ignore them. 8 For a discussion of a tower of effective theories relevant to the 
uncolored sector, see Ref. [21]. 

There are several things we have left out of Theory 1 because they do not affect color 
dynamics at leading order. One is a collision term in the Boltzmann equation (8a) due to 
hard collisions. The cross-section for such collisions is 0(a 2 ) and is overwhelmed, for color 
dynamics, by softer collisions that will be discussed in the context of Theory 2 below. 9 We 



7 Technically, W is the adjoint representation piece of the density matrix describing the color charges of the 
hard excitations, summed over the various species of excitations and integrated over the energy of excitations 
(for a fixed direction of motion v) . It is normalized in a way that simplifies the resulting equation. 

8 There are also other colored sectors besides the adjoint one represented by W. For each species of hard 
excitations, there are sectors corresponding to every irreducible representation contained in R x R, where 
R is the color representation of the given excitation. See Ref. [5] for a discussion. However, we do not know 
of any interesting physics that specifically couples to the other colored sectors except at sub-leading order 
in g. 

9 Hard collisions are relevant at leading order in g to physics in the uncolored sector, however. See, for 
example, Ref. [22]. 
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have also left out the soft modes of non-gauge fields such as scalars or fermions. These do 
not affect the color conductivity or the final soft dynamics (Theory 3) at leading order in g, 
provided the effective thermal mass of any colored scalars is large enough [large compared 
to g 2 T] so that they decouple at the very soft scales of interest for Theory 3. 

Theory 2. The w C fc < m limit of Theory 1 is a small frequency version of a theory 
originally written down by Bodeker [2], and is discussed in detail in Ref. [12]. In this low 
frequency regime, all time derivative terms in Theory 1 become negligible, and have been 
dropped, except for the one implicit in E (which is —dA/dt in A Q = gauge). Gauss' Law 
is replaced by the constraint (9b) due to the effects of Debye screening for k <C m. (See Ref. 
[12] for details and interpretation.) 

The most physically significant change is the introduction of the collision term 5CW 
into the linearized Boltzmann equation (9a). This term represents the effects of 2 — > 2 
collisions of hard particles via the t-channel exchange of what are called semi-hard gauge 
bosons (u ^ k ^ m). The leading log calculation of the linearized collision operator 8C 
was first made by Bodeker [2] and is explained in the language of collisions in Ref. [5]. One 
finds, 

5Cf(v) = (5C(vy)f(v')) v/ , (11a) 



ic(vy)» 7 W 

where 



J»(v-vQ-l / W) 



(lib) 



7 M ~C a c*T In (12) 

is the leading log contribution to the hard gauge boson damping rate due to exchange of 
gauge bosons having spatial momentum q > //. In other words, in a Wilsonian view of 
renormalization, \x is the ultraviolet cut-off imposed on k to define effective theory 2, and 
7(/i) is the damping rate due to gauge field fluctuations which have been integrated out 
and are no longer present in the effective theory. (In practice, it will be more convenient to 
perform renormalization by standard subtraction methods, and ji will be the renormalization 
scale rather than an ultraviolet cut-off.) The normalization of the angular delta function 
S 2 (v) is given in table II. 

The collision term in Theory 2 causes dissipation in the effective theory and so an appro- 
priately matched thermal noise term must also be present. The £ in (9a) is Gaussian noise 
with correlation [2] 

«£(v,x,tK(v',x',t')» = ^5C(v,v')5(t-t')5 (3) (x-x'). (13) 

A path integral representation of this effective theory may be found in Ref. [12]. 

Theory 3. This is Bodeker's effective theory, discussed earlier, and is obtained at leading 
log order by dropping \i below 7 (the momentum scale below which Theory 3 becomes valid) 
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in Theory 2, and then solving the W equations for fc C 7 [which allows one to drop v ■ D 
compared to 5C in (9a)]. Details may be found in Refs. [2,5]. We emphasize that, for the 
purpose of expanding in inverse powers of logarithms, we formally consider the momentum 
scale k ~ g 2 T of ultimate interest to be parametrically small compared to 7 ~ g 2 T ln^ -1 ). 

The decoupling of massive fields from low energy physics is very familiar in Euclidean field 
theory, and forms the basis for traditional applications of effective field theory techniques. In 
the case at hand, it is possible to eliminate the field W (while maintaining a local description 
in the effective theory) because one is dropping below the scale 7 characterizing the decay 
time of color correlations. This harks back to an old suggestion by Lebedev and Smilga [23] 
that 7 might cut off infrared divergences for some quantities. (See also Ref. [24].) As we 
shall see, this indeed happens for the color conductivity at next-to-leading log order, in a 
sense that we will make precise. 

It has long been known that the leading-order hard-thermal-loop effective interactions 
between gauge fields with momentum k ~ m ~ gT are non-local in space and time. That 
is, formally eliminating W from the Theory 1 Eqs. (8) generates non-local interactions for 
A. It should be emphasized that it is the inclusion of collisions which allows one to recover 
locality of effective gauge interactions for k <C 7, as in Theory 3. In the non-Abelian theory, 
arbitrarily small angle collisions will cause the color correlations of hard particles to fall off 
rapidly for times and distances large compared to the relevant collision time I/7, permitting 
one to treat such correlations as local in an effective theory restricted to k <C 7. Leading- 
order calculations using hard-thermal-loop interactions miss this effect because they do not 
incorporate the effect of collisions on the propagation of hard color fluctuations. In such 
a collisionless approximation, the disturbance created by probing the hard particles at one 
time does not decay, and can have an arbitrarily long-lasting effect on fluctuations measured 
at a later time — hence, non-locality. With collisions included, this non-locality is cut-off at 
the scale I/7. 



IV. SCALING AND SUBLEADING CORRECTIONS 

A. Review of static case 

One of the great advantages to using effective theories is that, besides simplifying the 
description of the relevant physics at a given scale, they provide a very clean mechanism for 
organizing, analyzing, and computing sub-leading corrections. For hot gauge theories, the 
quintessential example is the analysis of static equilibrium properties by Braaten and Nieto 
[25] using a sequence of Euclidean effective theories. (See also Ref. [26].) The important 
scales for static equilibrium physics are T, gT, and g 2 T; the scale 7, which will be crucial 
to our discussion of dynamics, does not appear. Below T, the effective theory is a three- 
dimensional Euclidean gauge theory coupled to an adjoint scalar A$ with mass m = O(gT), 
this mass being the manifestation of Debye screening. Below the scale m, the Aq field 
decouples, and one is left with a three-dimensional theory of unscreened magnetic physics: 



'eff 



: j'^'y-Vu ( i4 ) 



The term shown explicitly is simply renormalizable (in fact finite) 3-dimensional pure gauge 
theory. The dots represent an infinite sequence of interactions with higher and higher scaling 
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dimension (F 3 , F A , etc). These terms are infrared irrelevant, in the sense of the renormaliza- 
tion group, and their importance is suppressed by higher and higher powers of the scale k of 
the physics of interest, compared to the scale A = m ~ gT where the effective theory breaks 
down. So, for the study of non-perturbative physics, where k ~ g 2 T, terms in S c s with 
higher and higher scaling dimension correspond to corrections whose effects are suppressed 
by higher and higher powers of the ratio k/m ~ g. 

It is important to note that the coefficients of specific higher order terms may also contain 
explicit factors of g, which will cause their effects to be even more suppressed than indicated 
by their scaling dimension alone. 

Standard power counting of 3-dimensional gauge theory (14) gives the ultraviolet scaling 
dimension of the field A to be [A] = [rr] -1 / 2 (or equivalently [A:] 1 / 2 ). One typically rescales 
the field A — > T _1 / 2 A so that the coefficient of the kinetic term in (14) is dimensionless and 
so that counting scaling dimension is then the same as counting engineering dimension. In 
this paper, however, we shall avoid such redefinitions in order not to obscure the relationship 
between the gauge fields appearing in successive effective theories. 

Braaten and Nieto showed how, in principle, one can calculate the coefficients of the 
higher- dimensional terms in the effective theory (14), to whatever order desired, by carrying 
out perturbative matching calculations. Perturbative calculations suffice because the effec- 
tive theory is valid for scales k <C m, and because the physics of gauge field fluctuations is 
perturbative at the upper end g 2 T < C m of that range of validity. Roughly speaking, 
the idea is to pick a set of physical quantities that probe physics at a scale < A where the 
effective theory is valid, then to calculate those quantities in both the effective theory and 
the underlying theory, and then fix the (finite) set of unknown coefficients (at a given order) 
by requiring the answers to be the same. A refinement is required to make the calculation 
tractable, because long-distance physics in the effective theory (14) [and so also in the un- 
derlying theory] is non-perturbative! Conceptually, one imagines restricting the theory to a 
large but not too large box — one of size L which is large enough (L ^> A -1 ) for the effective 
theory expansion to be valid, but small enough (L <C (g 2 T)~ 1 ) that all the physics in the 
box is perturbative. The matching calculations can then be performed perturbatively. 

The "large box" temporarily introduced for the sake of matching can be replaced by any 
suitable infrared regularization, and in practice it is most convenient to use dimensional reg- 
ularization. Dimensional regularization is also used for ultraviolet regularization, if required 
to define the theory (at a given order) in the first place. 

B. The dynamic case 

Given that our goal is to expand in inverse powers of logarithms ln(g _1 ), while always 
working to leading order in powers of g, trying to refine our first two effective theories (8) 
and (9) by introducing new infrared-irrelevant interactions is unnecessary. We need to be 
able to use Theory 1 at the scale m in order to match onto Theory 2, but then irrelevant 
interactions are suppressed by powers of m/T ~ g, not by mere powers of inverse logs. 
We need to be able to use Theory 2 at the scale 7 in order to match onto Theory 3, but 
then irrelevant interactions are suppressed by powers of 7/m ~ g ln^ -1 ), which can again 
be ignored. The crux of an expansion in inverse logs, as far as the introduction of higher- 
dimensional interactions is concerned, is in the application of Theory 3 (Bodeker's effective 
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theory) to non-perturbative physics at a scale k ~ g 2 T <C 7. The influence of higher 
dimension irrelevant terms will then only be suppressed by factors of fc/7 ~ l/ln(g~ r ). 

Fortunately, a delightful simplification occurs, to be explained momentarily: effects of 
adding any higher-dimension corrections to Bodeker's effective theory (10) are suppressed 
by more than one power of the inverse logarithm. This means that Bodeker's effective 
theory is already perfectly adequate at next-to-leading log order, provided one calculates its 
one parameter (namely a) to next-to-leading log order. Previous numerical simulations of 
Bodeker's theory [3] can therefore be instantly extended to NLLO accuracy simply by using 
the improved value of a. 

Here is one version of the power counting analysis which demonstrates the required 
suppression of higher- dimension operators. It is possible to turn the Langevin equation 
(10) that defines the effective theory into a supersymmetric path integral form that allows 
standard renormalization analysis [12], which we review in the appendix. But one can get to 
the same result in a simpler, more cavalier way. The dynamical effective theory has to (and 
does) produce the same equilibrium distribution that is more traditionally analyzed with 
the Euclidean effective theories discussed earlier. We can therefore borrow the well-known 
result that the field A scales with distance as [A] = [x]^ 1 ^ 2 in the ultraviolet. That means 
that, at short distance, D — > V and Bodeker's effective theory becomes the free theory 
aE = V x B + C- This is just —a dA/dt = V x VA + £ in A = gauge, and implies that 
time scales as [t] = [x] 2 . This, of course, is just a re-iteration of (3). The lowest dimension 
terms that could possibly be added to the equation (10) for the effective theory, consistent 
with parity and gauge- invariance, are E x B and B x D x B. These terms have scaling 
dimension [x]~ 4 as opposed to those in Bodeker's equation (10), which are [re] -5 / 2 . The 
difference of 3/2 in scaling dimension means their effect is suppressed 10 by (fc/7) 3 / 2 . So, 
based simply on scaling dimension alone, the effect of possible higher dimension operators, 
when k = 0(g 2 T), is suppressed by at least [ln((? _1 )] _3 / 2 . This substantiates our claim 
that the form of Bodeker's equation (10) remains unchanged to next-to-leading log order. 
The effects of higher dimension terms may, of course, be even smaller that this estimate 
if there is explicit suppression in their coefficients. In fact, as discussed in the appendix, 
the coefficients of these operators necessarily include a factor of g which implies a further 
suppression of (g 2 T/'j) 1 / 2 . Hence, at the g 2 T scale, the actual relative correction due to 
higher-dimension operators is at least [ln((? _1 )] _2 . 



C. Implications for the meaning of a 

The above discussion allows us to provide an unambiguous definition of what we mean 
by the color conductivity at NLLO: we simply define it to be the coefficient of Bodeker's 
effective theory (10) at that order. This is a non-trivial issue because the conventional 
definition of conductivity for an Abelian gauge theory, in terms of the linear response of the 



10 It should be noted that this power counting argument gives the correct suppression factor when comparing 
terms that are either infrared irrelevant or marginal, and which do not have unnaturally small coefficients 
at the scale of the cutoff on the effective theory. It is the derivative terms in the uncorrected effective theory 
which define marginal scaling. 
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current j to a static, homogeneous electric field E, does not generalize in a meaningful way to 
a non-Abelian theory. This is because, in a non-Abelian theory, there is no gauge-invariant 
way to define the induced current j, nor is there any straightforward way to define the 
imposition of a "homogeneous" background electric field on top of the thermally fluctuating 
gauge field. 11 In contrast, Bodeker's effective theory is a local theory, and Eq. (10) is nicely 
gauge covariant with a a gauge-invariant numerical constant. 12 

Beyond NLLO, however, the definition of a as the parameter of the effective theory pre- 
sumably becomes ambiguous. That's because at next-to-next-to-leading log order (NNLLO), 
we need to introduce irrelevant infrared interactions into the theory, namely E x B and 
B x D x B, as discussed earlier. Such interactions render the theory non-renormalizable and 
so require the introduction of UV regularization to what was previously a UV-finite effective 
theory. That is entirely normal for an effective theory, but it means that the parameters of 
the theory will now be convention dependent, depending on the arbitrary choices of renor- 
malization scheme and renormalization scale. So, barring a conspiracy, we only have a clean, 
unambiguous definition of a at NLLO. It is entirely possible that color conductivity is only 
an approximate concept and that there is no natural, unambiguous definition beyond this 
order. 

An unambiguously defined conductivity is, of course, irrelevant to the success of the 
effective theory at any order, just as (the lack of) a convention-independent coupling constant 
is irrelevant to the applicability of familiar zero-temperature perturbation theory. 

V. THE RESULT FOR a 

The details of our calculation of NLLO color conductivity are given in Ref. [13], but we 
will summarize the main results here. Throughout, we use dimensional regularization in 
d = 3 — e spatial dimensions with gauge coupling /W 2 g. 

A. Matching Theory 1 to Theory 2 

Matching Theory 1 to Theory 2 involves calculating, to the requisite accuracy, the lin- 
earized collision operator 5C which results from the effects of gauge field fluctuations in 
Theory 1, in the presence of a perturbative infrared cut-off. Conceptually, one must then 
repeat this calculation in Theory 2 taking into account both the bare collision operator in 
Theory 2, and the residual effects of gauge field fluctuations which remain in Theory 2. A 
major virtue of using dimensional continuation for both infrared and ultraviolet regulariza- 



11 Actually, one can give a sensible meaning to a uniform non-Abelian gauge field [27]. But this does not 
generalize in a natural way to defining the addition of a uniform external field on top the thermal fluctuations 
in the gauge field. 

12 The alert reader may object that we earlier expressed uncertainty as to whether Eq. (10) is correct (and 
unambiguous) in gauges other than flow gauges. However, the Langevin equation (10) can be reformulated 
as a manifestly gauge- invariant path integral, as discussed in Ref. [12]. 
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tion is that, with this choice of regulator, the total collision operator in Theory 2 equals the 
bare operator; the residual gauge field fluctuations in Theory 2 have no net effect. 

The computation in Theory 1 may be carried out using any of the approaches that have 
been used for leading log results [2,5,3,16]. The kernel of the collision operator is given by 



6C(y, v') 



C A m 2 T 



2^V 



/ |M(v,v',q)| 2 \ ^(v-v')- / |M(v,v',q) 



(15) 



with the scattering amplitude M. which results from t-channel gauge boson exchange with 
hard-thermal-loop (HTL) self-energies included, 



A m 9fl J (27T)*- w " U 2 + n T (Q) Q 2 + n 



2tt5(Q • v) 2n5(Q ■ v') . 

(16) 



(P T and P L are transverse and longitudinal projection operators, respectively. See Ref. [13] 
for details.) The resulting form of the kernel <5C(v,v'), accurate to NLLO, is much more 
complicated that the leading-log expression (lib). Fortunately, the only portion of 8C which 
is actually needed for the determination of the NLLO color conductivity is the single matrix 
element 



7l = (tfSCv* 



(17) 



At leading- log order, this is the same as the hard gauge boson damping rate (6), but at 
NLLO these quantities differ. We find [13] that to leading order in g and all orders in 

i/inGr 1 ), 



7i(/i) = C A aT 



1 , / m 
- + ln — 



1 , 7r 
+ £ In - + ai 



(18) 



where we have written the result in terms of the MS scale /x = \^Arre ( 7E / 2 ^/i . The constant 
ai is given by the integral 



1 f 

ai = - y d\ I (l-A 2 ) 



1 argp T (A) argp L (A) Rg In /4(A) - lnp L (A) 



2 Imp T (A) Imp L (A) 



p^(A)-p L (A) 



1 

IAJ 



with 



PL(A) = l~ln| 



Pt(A) = - 



1 + A 



in 



A/ + T A ' 



Pl(A) 



.l-A 2 

(A more explicit integral representation for a± may be found in [13].) Numerically, 

ai = 0.323833 



(19) 

(20a) 
(20b) 

(21) 



As noted above, the corresponding calculation in Theory 2, when dimensionally regu- 
larized, is trivial, and matching Theory 2 to Theory 1 merely means setting the collision 
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operator equal to the result of (15) and (16), with the same matrix element 71 appearing in 
(18). This is the bare value of 71 in Theory 2; to obtain the renormalized value in minimal 
subtraction, simply drop the 1/e pole. 

It should be emphasized that our argument that the color conductivity a can be un- 
ambiguously defined at NLLO as the parameter of an effective theory does not apply to 
71 . Because Theory 2 is not ultraviolet finite, its bare parameters depend explicitly on the 
renormalization scale p, as seen explicitly in the above result for 71 (/x). This is completely 
analogous to, for example, the minimal-subtraction definition of quark mass in QCD. But 
a precise, convention independent definition of 7! is irrelevant (and unnecessary) for the 
correct evaluation of physical quantities. 13 



B. Matching Theory 2 to Theory 3 



Matching Theory 2 to Bodeker's effective theory at NLLO requires a one-loop calculation 
in both these theories. In Ref. [13], we describe how appropriate ratios of large time-like 
Wilson loops may be chosen as the physical quantities whose expectations are matched 
in order to determine the parameter a of Bodeker's effective theory. In practice, this is 
completely equivalent to matching the uj <C k, k — > behavior of the A self-energy IT 00 (cj, k) 
in Coulomb gauge. We find, 



1 - 



CpfiiT 

7 



+0 



.^(1/9)) 



(22) 



where p is again the MS scale, and the constant b is given by 



2 r°° , 
- / dp 

IT JO 



1 



So(p) - 



+ 



3p- 



■(l-E, 



4 v ( \ n P 

Is 2(p) ~ WTT) 



(23) 



Here, p represents the dimensionless magnitude |p|/T of a loop momentum, and £ m (p) are 
the dimensionless functions defined in bra-ket notation as 



E m (p) = (mm I (iv z p + 8C/i) 1 \mm), 



(24) 



where 5C is the linear operator defined by (11), the inverse is a non-trivial operator inverse 
(because 5C and v z do not commute), and \lm) represents the function Yi m (y). We give 
formulas for evaluating these functions numerically in Ref. [13]. The numerical result for b 
is 



2.8380 



(25) 



1 The matrix element 71 differs by a finite regularization independent amount from the hard gauge boson 
damping rate [13,24]. Many attempts have been made to formulate a gauge-invariant "pole mass" definition 
of the damping rate for non-Abelian theories, but these have so far led to a mire of confusion [29]. For 
discussion of the related (but quite distinct) issue in QED, see Ref. [30]. 
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C. Final result 



Combining Eqs. (18) and (22) yields our NLLO result for the color conductivity. The 
structure is clearest if we write an expansion for a -1 (the "color resistivity") rather than a 
directly. One finds 



_ x 3Ca«T 
a = 




+ C + 0\ 



>(l/<7). 



(26a) 



with 



C = |ln^ + a! + b = 3.0410- •• . (26b) 



Note that the 1/e divergences have canceled, as they must. Inside the logarithm of (26), 
7(/i) is to be understood as simply the leading-log formula 

7 (^) ~ C^T In (27) 

and \x should be chosen so that it is of order 7. One may easily verify that the \i dependence 
in the NLLO result (26) only affects that answer at order [ln(m/7)] _1 ~ [ln(l/(?)] _1 , which 
is beyond the order of this calculation. 

In order to use the NLLO result (26) in any practical calculation one must, of course, 
choose some particular value of \x [and ignore the unknown 0(1/ lng" 1 ) corrections]. In the 
absence of a full next-to-next-to-leading-log analysis, there is no clearly preferred procedure 
for determining an "optimal" value of /i. However, one reasonably natural choice is the 
"fastest apparent convergence" (FAC) scheme 14 which is to choose ji so that the next-to- 
leading order correction vanishes. In the present context, this amounts to choosing the scale 
/i FAC which satisfies 

/Ufac = e" C 7(/i FAC ) • (28) 
For this choice, the NLLO conductivity is simply 

<^o = — W— • (29) 



VI. CONCLUSION 

We have shown that at next-to-leading log order, no modification is required to Bodeker's 
effective theory for non-perturbative color dynamics other than the insertion of the correct 
NLLO value of the color conductivity. Using our NLLO result (26), one may instantly 



14 At our ncxt-to-leading-log order, one cannot use any minimal sensitivity criterion, since the NLLO result 
(26) has no stationary point in fj,. 
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generalize Moore's numerical result (7) for the topological transition rate of hot electroweak 
theory to NLLO: 



T « (10.8 ±0.7) 



in 



a T 



In 



m 



+ 3.041 + 



ln(Vs), 



(30) 



Moore, and others, have also obtained numerical results for the topological transition 
rate by using a more microscopic theory (analogous to a lattice version of what we called 
Theory 1) [31,32]. Moore also attempted to estimate the size of the NLLO correction to T 
by fitting the results of these simulations to the functional form 



r = k — 



9l 

rn 



<T 4 



In 



m 



+ 5 



(31) 



with the value of k fixed to 10.8, as determined from simulations of Bodeker's effective 
theory. This led to 5 3.6, with perhaps 20% uncertainty due to systematic errors [32]. 15 
This implies an estimate of 4.4 for the value of the \n(m/g 2 T)+5 factor in (31), which 
may be compared with the square bracket appearing in (30). This comparison is presented 
graphically in Fig. 1. The solid line is a plot of ln[m/'~f(fj,)]+C as a function of 7//X, for the 
specific case of electroweak theory with a single Higgs doublet and g 2 = 0.4. The dashed 
line indicates the value of 4.4 estimated in [32]. The arrow on the abscissa indicates the 
FAC point, where i(n)/n = e c ~ 20.93. 



In 
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FIG. 1. The value of the ln[m/7(/i)]+C factor, appearing in the inverse color conductivity (26) and 
in the topological transition rate (30), plotted as a function of 7//U, for electroweak theory with a single 
Higgs doublet and g 2 = 0.4. The dashed line indicates the value of 4.4 for which the NLLO result for the 
topological transition rate (30) (with the unknown yet-higher-order l/lng -1 terms neglected) agrees with 
independently determined results from more microscopic numerical simulations [31,32]. This value should be 
regarded as having a significant systematic uncertainty of perhaps 20%. The arrow on the abscissa indicates 
the FAC point where 7(/u)//x = e c ~ 20.93. 



15 Lattice artifacts exist in these more microscopic simulations (due, in particular, to the lattice dispersion 
relation allowing unwanted Cherenkov radiation), which cause them not to reproduce, precisely, the dynamics 
responsible for NLLO corrections to T. These effects have only been crudely estimated, but are a major 
part of the uncertainty in the estimate of 5. 
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The similarity between our NLLO result and the value inferred from numerical simula- 
tions is remarkable. There was no obvious reason a priori why it should be a reasonable 
approximation to treat logarithms of the gauge coupling [that is, ln(#/g)], as large for 
physical values of the coupling. However, one sees from Fig. 1 that even the most naive pre- 
scription of picking \x exactly equal to the scale 7 produces a result within 20% of Moore's 
numerical estimate (which itself has an uncertainty estimated at 20%). The naive prescrip- 
tion of letting n vary by an entire order of magnitude about 7 (from O.I7 to IO7), and 
taking the variation as a crude estimate of the uncertainty from leaving out higher-order 
corrections, still keeps the result within roughly 35%. The FAC value of ^(/i) / /j, suggests 
that perhaps one should pay more attention to the range /i < 7, where the variation with \x 
is even smaller and the agreement with Moore's numerical estimate even better. 

The close agreement between the FAC value of [i and the precise point where the curves 
of Fig. 1 cross is striking, but probably fortuitous given the uncertainty in the numerical 
simulation value. Nevertheless, perhaps the characteristic size of neglected corrections really 
is e~ C 7, and not just 7 (as one might naively expect). If true, this would mean that (for 
electroweak theory) the expansion in inverse logs actually has a respectably small expansion 
parameter of about 0.25. To confirm (or refute) this expectation, one would need to cal- 
culate to yet higher order in the expansion in inverse logarithms. As discussed earlier, this 
presumably requires adding new operators to Bodeker's effective theory, and this in turn im- 
plies that new numerical simulations would be required. Such additional higher-dimension 
operators will also cause the effective theory to no longer be UV finite, implying that a 
careful matching calculation for the lattice-regularized theory would now be required. 
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APPENDIX A: POWER COUNTING WITH FUNCTIONAL INTEGRALS 



In the main text, we discussed how to count the relative scaling dimensions of operators 
appearing in the Langevin equation. We swept under the rug issues of how to do the same 
for terms involving the noise. To be systematic, it is much more convenient to recast the 
problem in terms of a functional integral, where dimension counting of operators is more 
familiar. Also, in the functional integral version, noise has been integrated over, and need 
not be discussed separately. 

The functional integral representation of Bodeker's theory is well known as the functional 
integral representing stochastically quantized 3-dimensional non-Abelian gauge theory. See, 
for example, Ref. [33] for a review, or [12] for a discussion in the present context. In A = 
gauge, one form is 



Z = J[DA(x,t)] J[A] e- s ^ 



with 



dt (Fx 



dt 



(Ala) 



(Alb) 



and where J [A] is a Jacobian factor that equals one in dimensional regularization. It is 
useful to rewrite the action in terms of the magnetic energy density, 



V(A) 



±B a • B a , 



(A2) 



as 



S = 



AaT 



dt d 3 x 



dA dV 

a 1 

dt dA 



(A3) 



There is a hidden supersymmetry of this functional integral which restricts the form of 
interactions which can be added to it. The functional integral can be rewritten in superfield 
notation as 



with 



1 

T 



where <fr is the superfield 



Z = 



= J im 



-S[*] 



d9 d6 dt d 3 x 



aV t <f>-V t <$> + V(<!>) 



(A4) 



(A5) 



$ = a + 9c + c9 + xee. 



(A6) 



The supersymmetry only involves time, rather than space-time, and the SUSY time deriva- 
tives are 



f> t = de, 



V t = de - Odt . 



(A7) 
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In (A6), A is an auxiliary field, and c and c are ghost fields, the integral over which produces 
the Jacobian factor J [A] mentioned earlier. 

The scaling dimensions of all the fields may now easily be read off from the supersymmet- 
ric action: [t] = [x] 2 , [9] = [9] = [x], and [<&] = [re] -1 / 2 . (Keep in mind that for Grassmann 
integration, d9 has the inverse dimension of 9, and so [d9] = [dO] = [x]^ 1 .) One may now 
analyze the dimensions of what possible irrelevant interactions can be added to the super- 
symmetric action (A5) consistent with gauge invariance, parity, and supersymmetry. 16 This 
provides a more familiar way to do the power counting analysis than that of the main text. 
For instance, consider the interaction BD 2 B, where B and D are to be understood as the 
normal non-SUSY expressions for the magnetic field and the gauge-covariant derivative, but 
with A replaced by the superfield The interaction BD 2 B has scaling dimension [x]~ 5 , 
as opposed to the other terms in the supersymmetric Lagrangian, which have the marginal 
scaling [a;] -3 . So the effects of this term should be suppressed by at least (fc/7) 2 , since 7 
serves as the cut-off scale A for this theory. Another example is V t & ■ V t <& x B (still keeping 
to A = gauge), which has scaling dimension [a;] -9 / 2 , and so is suppressed by (fc/7) 3 / 2 . In 
fact, the coefficient of this operator must also contain a factor of g, because the underlying 
theory is unchanged if g — > —g and A — > —A (or <fr — > — <&). The effective theory must also 
respect this invariance. However a suitable rescaling of variables in Theory 2 (see Ref. [13]) 
shows that the theory only depends on the scale 7 plus the dimensionless combination g 2 T 
Hence, a factor of g in the coefficient of an induced higher dimension operator implies that 
the effects of this operator will have an additional suppression by (g 2 T/'j) 1 ^ 2 . Therefore, the 
effects of the interaction T>t*& • T>t& x B must be suppressed by (g 2 T/'j) 1 ^ 2 (fc/7) 3 / 2 . 

A little thought shows that the net result is that the only interactions which can be 
added to the action (A5), consistent with its symmetries, are either (i) irrelevant operators 
whose effects, at the scale g 2 T are suppressed by at least two powers of g 2 T/^, or (ii) total 
derivatives [such as (V t &) • D x B = D t V(<&)] which have no effect on the dynamics. 



16 Although it is not required for our discussion, it is worth noting that time-reversal invariance also con- 
strains the possible irrelevant interactions which may appear in the effective theory. This may be surprising 
at first sight, since the Langevin equation (10) defining the leading-order effective theory is dissipative and 
manifestly violates time reversal invariance. Nevertheless, this theory generates time-dependent equilibrium 
correlation functions which do respect time-reversal symmetry. Formally, this is easiest to see from the 
functional integral (Al). The action density appearing in (Alb) is not invariant under time reversal. How- 
ever, after multiplying out the square, it is only the cross term which is time reversal odd. And this cross 
term is a total time derivative, J d 3 x a (dA/dt) • D x B = (d/dt) J d 3 x crV(A). Hence, the action (Al) is 
in fact time-reversal invariant up to boundary terms at t = ±00, which are irrelevant as far as equilibrium 
properties (including time-dependent correlation functions) are concerned. This invariance (up to boundary 
terms) must remain true in the presence of higher-dimension irrelevant interactions. For related discussion 
in terms of the Fokker-Planck equation, see Ref. [34]. 
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